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Abstract 

We propose a scheme for generating vibrational pair coherent states of the 
motion of an ion in a two-dimensional trap. In our scheme, the trapped ion is 
excited bichromatically by three laser beams along different directions in the 
X-Y plane of the ion trap. We show that if the initial vibrational state is given 
by a two-mode Fock state, the final steady state, indicated by the extinction 
of the fluorescence emitted by the ion, is a pure state. The motional state of 
the ion in the equilibrium realizes that of the highly-correlated pair coherent 
state. 
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A variety of generalized coherent states has been constructed to describe different physical 
phenomena Mathematically, the constructions of these sets of generalized coherent 
states are associated with particular Lie groups. The Glauber coherent states, defined as 
the right eigenstates of a single-mode boson annihilation operator are associated with the 
Heisenberg group. Beyond these canonical coherent states, particular generalized coherent 
states associated with the noncompact SU(1,1) group have been extensively studied 
According to Barut and Girardello these associated coherent states are defined as the 
eigenstates of the SU(1,1) lowering operator, whereas in the definition given by Perelomov 
10], they are generated by the SU(1,1) analogue of the displacement operator. These two 
sets of coherent states are different though they are closely related to the SU(1,1) group. 

Regarding the two-mode boson realization of the SU(1,1) group, a special set of coherent 
states of the Barut-Girardello type known as pair coherent states (PCS) can be formulated 
If a (a^) and b (b"^) denote two independent boson annihilation (creation) operators, then 
ab (a^6^) stands for the pair annihilation (creation) operator for the two modes. The pair 
coherent states |^, q)pcs are defined as eigenstates of both the pair annihilation operator ab 
and the number difference operator Q = a^d — 6^6, i.e.. 



where ^ is a complex number and q is the "charge" parameter which is a fixed integer. 
Furthermore, the PCS can be expanded as a superposition of the two-mode Fock states. 



where J\fq = [\^\ ''/q(2|^|)] is the normalization constant (Jg is the modified Bessel 
function of the first kind of order q). Pair coherent states were introduced by Horn and Silver 



which possess prominent nonclassical properties such as sub-Poissonian statistics, correlation 
in the number fluctuations, squeezing, and violations of Cauchy-Schwartz inequalities [§. 
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abl^, q)pcs = CIC, g)pcs , Q\^, q)pcs = q\^, q)pcs , 



(1) 




(2) 



^ to describe the production of pions and applied to other problems in quantum field theory 
In quantum optics, PCS are regarded as an important type of correlated two- mode state. 



The experimental realization of such nonclassical states is of practical importance. Agar- 
wal suggested that the optical PCS can be generated via the competition of 4-wave mixing 
and two-photon absorption in a nonlinear medium. This is the only scheme proposed to gen- 
erate PCS known to us. In this Rapid Communication we propose a scheme to generate 
vibrational PCS of a trapped ion. Recently, due to remarkable advances in laser cooling of 
a trapped ion it has become possible to realize nonclassical states of the centre-of-mass 
(CM) motion of a single trapped ion. An ion confined in an electromagnetic trap can be 
regarded as a particle with quantized CM motion moving in a harmonic potential. Exciting 
or deexciting the internal atomic states of the trapped ion by a classical laser driving field 
changes the external states of the ion motion, as atomic stimulated absorption and emission 
processes are always accompanied by momentum exchange of the laser field with the ion. 
If both the vibrational amplitude of the ion is much smaller than the laser wavelength, i.e., 
in the Lamb-Dicke limit [0], and the driving field is tuned to one of the vibrational side- 
bands of the atomic transition, then this model can be simplified to a form similar to the 
Jaynes-Cummings model (JCM) |^ in which the quantized radiation field is replaced by the 
quantized CM motion of the ion. As the coupling between the vibrational modes and the 
external environment is extremely weak, dissipative effects which are inevitable from cavity 
damping in the optical regime, can be significantly suppressed for the ion motion. This 
unique feature thus makes it possible to realize cavity QED experiments without using an 
optical cavity. Following this approach, nonclassical vibrational states of the trapped ions 
such as Fock 0, squeezed and Schrodinger cat states |irT| , p!2| have been proposed and 



observed |TI 

Consider the quantized motion of a two-level ion of mass M which is trapped in a two- 
dimensional (2-D) isotropic harmonic potential characterized by the trap frequency v. The 
creation (annihilation) of vibrational quanta in the X and Y directions is described by the 
operators a (a^) and h (6^) respectively. The position operators are given by a; = A(a+a^)and 



y = A(6 + 6^), where A = y ^/2z/M. In our scheme, which we sketch in Fig.|T], the ion is driven 
bichromatically by three laser beams in the X-Y plane. The first two lasers are both tuned 



to the second lower vibrational sideband and applied to the ion along directions with an 
angle 7r/4 and 37r/4 relative to the X-axis, respectively. The third laser which drives the 
ion along the X-axis is resonant with the atomic transition frequency. In the rotating-wave 
approximation, the Hamiltonian describing the coherent evolution is 



H = hiy(^ 



a' a 



hujQ 



VE^-^ {x,y,t)a_ + h.c. 



(3) 



where the first two terms describe the free evolution of the external and internal degrees of 
freedom of the ion and the last indicates the atom-field interaction. The operators o"+ and 
a- are raising and lowering operators for the two-level ion obeying commutation relations 
[o"_|_,(T_] = and [o"2,(3"-|-] = ±2d'j-. The transition in the two-level ion is characterized by 
the dipole matrix element V and the transition frequency uq. The negative frequency part 
of the classical electric driving field is given by 



_j_^^g«('^o*-fco£+</'o) 



(4) 



where 2Ej and 0j indicate the amplitudes and phases of the driving lasers respectively. We 
have introduced new position operators x' and y' which are related to x and y by a 7r/4 
rotation in the X-Y plane, 

' 1 i\(A 



1 
V2 



V 



-1 1 



(5) 



so that for the creation (annihilation) operators A (A^) and B (B^^) , defined in the X' and 
Y' directions respectively, we obtain the following transformation 



1 

71 



( 1 iUa\ 



-1 1 



(6) 



According to Vogel and de Matos Filho |TT| , p!^ , if the ion is in the resolved sideband limit 
and the driving laser is resonant with one of the vibrational sidebands, then the ion-laser 
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interaction can be described as a nonlinear JCM. Following Refs. [^|T4|, the Hamiltonian 
of Eq.(3) can be written in the interaction picture as 



,eoj!(j + 2)! 



(T_ + h.c. , 



(7) 



where Qj = —VEj/h are the Rabi frequencies, and the Lamb-Dicke parameter rj = kX has 
been defined assuming ^2 — = ^• 

As the damping of vibrational quanta can be significantly suppressed in an ion trap, the 
dominant decay process is the spontaneous emission from the two-level ion, and the time 
evolution of the system in the interaction picture can be described by a density operator p 
obeying the master equation HJTT[ 



dp 
'dt 



Hh p\ + 2 (20-- ^0-+ - o-+o-_/3 - o-_6-+p) 



(8) 



where T is the spontaneous decay rate of the excited state of the ion, and the modified 
density operator 

1 



' = AJ-^J-^ 



(9) 



accounts for the momentum transfer in the X-Y plane due to spontaneous emission, where 
W {u, v) describes the angular distribution of the spontaneous emission. In the Lamb-Dicke 
regime, ?7 -C 1, the master equation, Eq.(8), can be well approximated by the first order 
expansion in rj. In this case, only the leading terms, i.e., j = in Hi, are considered, and 
g is replaced by p. Setting ^1 = ^2 = ^, and 0i = 0, 02 = tt, the master equation Eq.(8) 
simplifies to 



dp 
'dt 



Hip] + 2 (2(T-/3(T+ - (3-+(T_p - (T_(T+/3) 



(10) 



with the effective Hamiltonian H{ given by 



H[ = a ah — ^ 



a+ + h.c. 



(11) 



where a = hflr]'^ exp {—r]'^ /2), and ^ = QqQ ^t] ^exp(— 20o)- Disregarding the constant 
driving term ^, the effective Hamiltonian of Eq.(ll) is identical to the non-degenerate two- 
mode JCM 

For a master equation of the form of Eq.(lO), the steady-state solution ps is a pure state 

Ps = |^)|^)(^|(^?|, (12) 

where \g) is the atomic ground state and defines the vibrational state of the ion. When 
the system reaches the steady state dps/dt = 0. Assuming a steady state of the form above 
gives [H{, ps] = 0, as the dissipative term on the right hand side of Eq.(lO) vanishes in 
this case. Consequently, a sufficient condition for us to generate the vibrational steady 
state lip) is to ensure ps satisfies the restriction [H{, ps] = 0, or equivalently dblip) = ^{ip) ■ 
Evidently, this is identical to our first definition of a PCS given above in Eq.(l). In order to 
generate a fully characterized PCS, however, we need another constraint to ensure that the 
final vibrational state remains an eigenstate of the number difference operator Q [Eq.(l)]. 
This can be achieved by properly choosing the initial vibrational state. In view of the fact 
that [db, Q] = 0, it follows that the number difference is a constant of motion in processes 
involving simultaneous pair annihilation or creation. Thus for an initial vibrational state 
\ipo), if the conserved "charge" condition is initially satisfied, i.e., Qlipo) = lli^o), then this 
condition holds for the whole time evolution of the system described by Eq.(lO). Regarding 
the feasibility of an experimental realization, it seems that the two-mode Fock state \iIjq) = 
\m+q, m)-p is best suited to fulfil the conserved "charge" condition. Without loss of generality 
we may set the initial vibrational state as iV'o) = lo', 0)f {q > 0) which, according to recent 
work, can be prepared with very high efficiency ||13|| . 

In order to identify the validity of our analytic argument and to gain insight into the 
transient behaviour we have solved the master equation, Eq.(lO), numerically employing a 
Monte-Carlo state- vector method . Our numerical analysis was performed using a high- 



order unravelling technique ||T^ in a finite (truncated) Fock state basis with a cut-off chosen 
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such that an increase of this cut-off does not alter the result of our integration. 

We find that the system evolves from an initial product of a two-mode Fock state and a 
superposition of the internal states of the ion into the pure steady state [Eq.(12)] as expected. 
Fig.0 depicts the excitation number distribution in the two vibrational modes at different 
times, showing the system as it evolves from the initial pure state |\I'o) = |e) ® |7, 6)f into its 
steady state. The steady-state excitation number distribution [Fig.H(d)] is indistinguishable 
from that of a PCS with ^ = 2 and q = 1- This is the expected result since the particular 
choice of the initial state determines the conserved "charge" q = (\l/o|Q|\l/o) = 1 , and from 
Eqs.(l) and (11) the second number defining the PCS is determined by the driving term ^ in 
the effective Hamiltonian, which in the example shown is ^ = 2 . We verify that the steady 
state is pure by calculating the trace of the square of the density operator Tr(pg) from the 
steady state obtained in our numerical integration. We find Tr(pg) = 0.9997. Fig.|^ shows 
the transient behaviour of the internal state of the ion. We depict the time evolution of the 
inversion, (a^), [Fig.|^(a)] and the imaginary part of the polarization, i (6'_ — a^) , [Fig.|^(b)]. 
We note the timescale for the system to reach the steady state is much longer than it would 
take the same ion to reach its steady state in a resonance fluorescence experiment. 

One may examine the existence of the PCS by observing the collapses and the revivals 
of the atomic inversion. This is illustrated by switching off the carrier field suddenly once 
the system has reached the steady state. In this case, the system is equivalent to the non- 
degenerate two-mode JCM interacting with a PCS when the ion is initially in its ground 



state. The time evolution of such a system has been investigated by several authors |T8|. In 
particular, when q = the Rabi oscillation of the atom exhibits exactly periodic behaviour. 
However, as the atomic decay is essential to the present treatment and cannot be ignored, 
the collapses and revivals of the Rabi oscillation of the trapped ion are substantially different 
from the previous results. One possible way to suppress the complications caused by the 
atomic decay is to increase the intensities of the driving lasers so that the ratio F/a is lowered. 
Thus, in the short time regime the influence of spontaneous emission can be eliminated. 
In conclusion, we have proposed a scheme for the realization of pair coherent states of 



the CM motion of a trapped ion. In our scheme, three laser beams, one of which is tuned 
to the carrier frequency and the other two to the second lower vibrational sideband, are 
used to drive the ion trapped in a 2-D isotropic harmonic potential well. In appropriate 
limits, the system will relax to a steady state due to the spontaneous emission from the 
ion. If the vibrational state of motion of the ion is initially prepared in a Fock state, then 
the steady state of the system is a pure state given by a product of the atomic ground 
state with a PCS of the vibrational motion. In this case, the two parameters, ^ and g, that 
characterize the PCS are determined by the intensities and phases of the driving lasers and 
by the number difference between the two vibrational modes determined by the initial Fock 
state respectively. 
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FIGURES 




FIG. 1. 

Configuration to generate vibrational PCS. The carrier field drives the ion on resonance 
along the X-axis. The other two lasers are both resonant with the second lower vibrational 
sideband. They drive the ion along the X' and y'-axis having phases (f)i — and 02 = tt 
respectively. 
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FIG. 2. 

We depict the excitation number distribution P(n,m) in the two vibrational modes at 
different times. The system evolves from (a), the initial pure state |^o) — \e) <Si |7, 6)f 
(n = 0.0), through (b) and (c), two intermediate states (Tt = 125.0) and (Tt = 500.0), into 
(d), its steady state which is the pure product of the internal ground state of the ion and a 
vibrational PCS: |*s) = \g) (g) |2, l)pcs {Tt = 2000.0). The data shown has been obtained 
from a Monte-Carlo simulation which included 1000 trajectories and was performed in a 
truncated Fock state basis (|0, 0)f...|20, 20)f). Parameters: a = 0.2, ^ = 2.0, F = 10.0 . 
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FIG. 3. 

Time evolution of the internal state of the two-level ion as the system evolves into its 
steady state, (a) shows the inversion, ((Tz), and (b) the imaginary part of the polarization, 
— (T_|_) . The real part of the polarization, {a. + 0"+), remains zero at all times. The two 
insets show the evolution in the time interval Tt = 0.0 to Tt = 20.0 . Parameters and data 
used are the same as in Fig.||. 
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